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We show that the (7(2) family of point interactions on a line can be utilized to provide the U(2) 
family of qubit operations for quantum information processing. Qubits are realized as localized 
states in either side of the point interaction which represents a controllable gate. The manipulation 
of qubits proceeds in a manner analogous to the operation of an abacus. 
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The idea of representing numbers as spatial locations 
is probably as old as mathematics itself. Even today, a 
functioning example of locational representation of num- 
bers is found in the form of a calculator known as abacus, 
in which digits, or bits, are stored as locations of one- 
dimensionally mobile objects. The purpose of this note 
is to present a quantum version of the abacus, a model 
system in which qubits are realized as spatially localized 
states in a one-dimensional space that is divided into two 
identical regions by infinitely thin barrier. The key to this 
realization is the concept of quantum point interactions 
HiHi, which is a mathematical abstraction of thin 
barriers acting on quantum particles. It is shown that the 
U(2) parameter family of all possible quantum point in- 
teractions, with the use of time evolution operators, can 
be turned into the U (2) structure of operations on a two 
dimensional qubit Hilbcrt space. Specifically, the S 2 sub- 
family formed by the scale-invariant point interactions is 
turned into the Bloch sphere of qubit transformations. 

The standard realization of a qubit Hilbert space em- 
ploys a system with two eigenstates, a spin one-half sys- 
tem being its archetype p]- Being purely quantal in na- 
ture, the system realized this way tends to demand very 
delicate handling for qubit operations that often becomes 
quickly impractical in many-qubit settings. In contrast, 
the quantum abacus model presented here possesses a 
semiclassical nature, as the system employs a qubit ba- 
sis comprising infinitely many eigenstates,which is made 
possible by utilizing quantum caustics |g,|2j- As a result, 
the system can acquire robustness and spatial extendabil- 
ity, which are highly desirable for quantum information 
processing beyond the level of laboratory experiments. 

The system we consider for qubit operations consists 
of a particle of mass m moving on a line (-co < x < oo) 
under a point interaction at x = and a background 
harmonic oscillator potential with frequency lu (FIG. 1). 
To describe a state ip(x) of the system, it is convenient to 
introduce ?/>+(ie) = ip{x) and ip-(x) = tJj(-x) for x > 
and thereby represent the wavefunction on the positive 
and negative sides separately, using the two-component 
wavefunction 



U 



FIG. 1: Our setting for qubit operations: A harmonic oscil- 
lator potential with a tunable point interaction at the center 
that provides a connection condition with a U(2) matrix U. 



The state of the particle is governed by the Hamiltonian 

(2) 

and is subject to the connection condition [sHiHTol 



2m dx 2 2 



d* 

{U-I)*(0) + iL (U + I) — {0) = 0, (3) 

which is specified by the matrix U € U(2) characterizing 
the point interaction. Here, / is the two-by-two unit ma- 
trix, Lq is a length constant required on dimensional 
grounds, and 'J and d^/dx are evaluated in the limit 
x +0. 

Among the point interactions described by J2J, the fa- 
miliar example of (5-interaction, which induces disconti- 
nuity in the derivative of the wavefunction, is only a spe- 
cial one-parameter subfamily. In fact, the U(2) variety 
comprises such exotic point interactions that cause dis- 
continuity in the wavefunction itself, as well as ones that 
admit a constant transmission probability for all particle 
energies. Explicit construction of these highly singular 
point interactions has been achieved in terms of sin gula r 
short-range limits of simple known interactions Il2| . 
To illustrate by some examples, we mention that the 



= 0, 



x>0. 



(1) 



identity matrix U — I results in 7^(0) = ^—(0) 
signifying a point interaction for an inpenetrable barrier 
at x = with Neumann boundaries at its both sides. 
Similarly, the negative identity matrix U = —I gives 
ip+(0) = "0-(O) = 0, another inpenetrable barrier with 
Dirichlet boundaries. Furthermore, it is the Pauli matrix 
U = o\ that provides the free case (i.e., no point inter- 
action) with the smooth connection conditions ^+(0) = 
-0_(O), 7^(0) = — -xjp- (0). Another important example 
is the case of Hadamard matrix U = (a\ + a^)/^/2. It 
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turns out that this point interaction admits the constant 
transmission probability one-half irrespective of the en- 
ergy of the incident particle @. 

A notable property of point interactions is that, given a 
state ^ that fulfils Q with U, we obtain, with an SU (2) 
matrix V, a new state VP' = V^f that satisfies © with 



u' = vuv~ 



(4) 



The transformations by various V provide a link between 
different point interactions, and lead to several salient 
properties among the family of point interactions 0, llf)j | . 
We utilize three of them, which are available when the 
background potential is parity invariant. 

The first is that if \t is an energy eigenfunction then 
V?' is also an energy eigenfunction with the same energy. 
This implies that the two systems with their point in- 
teractions related by the conjugation Q share the same 
energy spectrum. Since those U which fall into the same 
conjugacy class under Q form a sphere S 2 , all such 
spheres in U(2) are isospectral subfamilies. 

The second property concerns the time evolution under 
point interactions U, described by the transition operator 
U t {U) 



e~i Ht as 



(5) 



Suppose Ut(U) is the transition operator under the point 
interaction U , then states under U' evolves with the tran- 
sition operator 



u t (u') = vu t (u)v- 



(6) 



in correspondence to (@J. 

The third is the fact that the conjugation (QJ can be 
used to link two physically distinct point interactions, 
one separating (representing an inpenetrable barrier) and 
the other nonseparating, relating thus a barrier totally 
reflecting from both directions to a barrier which admits 
probability flow across x = 0. Indeed, since any U(2) 
matrix U can be decomposed in the form 



U = V~ l D V, 



(7) 



with V G SU(2) and a diagonal matrix D, the conjuga- 
tion (J3J with the same V amounts to the diagonalization 
of U into U' — D. If we parametrize D as 



D 











(8) 



with < 9± < 2ir, then the connection condition 10 for 
U = D reads 



V>±(0) + L cot%^(0) = 0, 

2 dir. 



(9) 



showing that the positive and negative sides are sep- 
arated physically. Note that these D in © form a 
U(l) x U(l) torus subfamily in U(2). Since the spec- 
trum is unchanged under conjugation, these D index the 



possible energy spectra of point interactions [ljj. Visu- 
ally, if the barrier represented by the point interaction is 
regarded as a gate, and if U is non-diagonal allowing for 
transmission, then the conjugation U — > D (or its inverse 
D — > U) is equivalent to opening (or closing) the gate. In 
general, the conjugation U — > U' by alters the trans- 
mission property of the gate while preserving the energy 
spectrum. Hence, the conjugation provides a very useful 
means to control the gate, and plays a vital role in the 
present proposal of quantum information processing. 

Our realization of qubit operations employs the class 
of point interactions given by 



U 



V e SU(2). 



(10) 



This conjugacy class possesses no finite scale parame- 
ter and provides (together with U = ±7) the subfamily 
of scale invariant point interactions in the U(2) family. 
Moreover, the class of the unitary matrices l|10(l leads 
to the Bloch sphere when implemented by time evolu- 
tion. This can be made more explicit by parametrizing 
V = e's^e''" 173 with angles < (i < tt and < v < 2tt 
to obtain U = <r(c) where a(c) is a vector in the linear 
space spanned by the Pauli matrices, 



(11) 



with (ci,C2,C3) = (sin [i cos v, sin (i sin v, cos /i) . In par- 
ticular, if fj, = 5, v = we have U = <j\ which yields 
the smooth connection condition. The system then be- 
comes the standard harmonic oscillator with the energy 
eigenvalues 



e n = (n+ l/2)Sw, n = 0,1,2,.. 
and the eigenfunctions 



u n (x) 

(-l) n U n (x) 



(12) 



(13) 



where u n (x) are the familiar harmonic oscillator wave- 
functions consisting of Hcrmite polynomials. Since the 
harmonic oscillator potential is parity invariant, all U of 
(fTUfl share the same energy spectrum by the second prop- 
erty of conjugation mentioned before. The conjugation 
also allows us to write the eigenstates for the point 
interaction 03 as 



CT 3 — e l 4 (72 ^) a ' 1 = 



V2u, 




-V2u n 



(n = even) 
(n = odd). 



(14) 



We now consider the time evolution of the system un- 
der U for a time step r, which we set to be the half period 
of harmonic oscillation, r = n/u). For the case of point 
interaction 03, we find 



U T (* 3 )$% = e"^<C = -ia s $%, (15) 
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FIG. 2: Our qubit basis of states |0) and These are 
provided by a state of arbitrary profile localized on one side 
of the barrier and the other state given by the mirror image 
of the former. 
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for which the relation e _ "K e " T = — i(— 1)™ is used. Since 
an arbitrary state 'F is expressible as a linear combination 
of , its evolution under the point interaction 03 in half 
period r is also given by U T {az)^ = —icr^. Thus we 
obtain a remarkable expression 

W r (cr 3 ) = -i 0-3, (16) 

namely, apart from the constant phase —i, the time evo- 
lution operator under the point interaction cr 3 is given by 
03 itself. With the conjugation @, the time evolution 
of a state in a system under the scale invariant family of 
point interaction U is given by the evolution operator 

U T {U) = -iU, (17) 

which is (up to —i) exactly the matrix that character- 
izes the point interaction. In brief, the spatial property 
encoded in a point interaction U manifests itself in the 
time evolution U T {U). 

Note that the present realization of unitary operations 
associated with the Bloch sphere (|1(J|) is available for an 
arbitrary state, not just for a particular pair of eigen- 
states. This allows one to consider a qubit space spanned 
by a state with an arbitrary profile localized on one side 
of the barrier and its mirror state. A convenient choice 
for the qubit basis |0) and |1) is 

*+ = ( fi o ] ) " 1°) > *"=(/(*)) "|1>, ^ 

where f(x) is an arbitrary real function having the nor- 
malization property J °° f(x) 2 dx = 1 (FIG. 2). Since the 
spatial profiles of the qubits |0) and |1) are mirror sym- 
metric, any mixtures of |0) and |1) that emerge as a result 
of the time evolution U T (U) have profiles that are analo- 
gous to each other apart form the different scaling factors 
on both sides of the barrier (FIG. 3). In other words, the 
state returns to the qubit space after the unitary evolu- 
tion associated with the Bloch sphere for a time step t. 

To implement an arbitary unitary operation beyond 
the Bloch sphere, we recall that a generic element U £ 
U(2) admits the decomposition iJBJ. It is easy to show 
that, instead of conjugation by V, a Bloch sphere element 
a(c) with properly chosen Ci may also be used to write 
U = cr(c) D a(c) [l(|. The diagonal element D, in turn, 
can be expressed as the product of two Bloch elements 
and an overall phase D = e l * a (a) cr(b). This is achieved 



FIG. 3: Our qubit operation is achieved via the time evo- 
lution of the initial state that has an arbitrary profile. For 
the time step r= tt/ui, the point interaction at the origin 
is switched from the closed gate U = -/ to an appropriate 
open gate U. The time evolution on the initial state, here 
illustrated by |0), yields the unitary operation U used for the 
opening gate by 1171 . The profile (modulo the scale) is kept 
after each time step. 

by using a set of vectors a(a) and a(b) of the form 
which arc orthogonal to each other in the plane perpen- 
dicular to the third axis, as can be confirmed from the 
elementary identity a(a) a(b) = (a ■ b)I + i J2i( a x ^)i a i- 
Since the phase e 1 ^ may be cancelled out by adding a 
constant potential to the Hamiltonian (J2J, we find that 
any unitary evolution U can be performed essentially by 
the successive application of corresponding time evolu- 
tion operators belonging to the Bloch sphere in up to 
four steps, each of which is realized by the time evolu- 
tion for the half period r. We conclude, therefore, that 
by tuning the point interactions appropriately, all the 
U(2) operations required for quantum processing can be 
implemented in the qubit space spanned by and 

There is a more direct method to implement the di- 
agonal operation D, if we do not stick to the operation 
solely by the point interactions. In order to obtain sepa- 
rate phases for |0) and |1), we can simply apply an extra 
constant potential on either side of the gate in addition 
to the overall constant potential, while keeping the gate 
closed by setting U = ±1. 

Let us now consider the trasition between the classical 
and the quantam regime in our abacus operation. The 
classical regime arises when we have a state of a sharply 
localized wave packet on which we consider only 'classi- 
cal' gates U = —I and U = cr\ (up to the sign). In each 
time step r, the wave packet, which can now be inter- 
preted as a classical bead, bounces back at the gate and 
stays in the original side under the closed gate U = —I, 
or travels to the other side with the open (or NOT) gate 
U = <7i. Obviously, this is a microscopic realization of 
the classical abacus. 

If, on the other hand, we bring all possible quantum 
gates U £ U(2) in use, bilocal wave packets inevitably 
come into play, as seen in FIG. 3. We then recover the 
full quantum abacus. Thus, if we can implement the 
full set of quantum gates on a larger scale which is 
possible as long as the harmonic potential is maintained 
accurately to the extent that admits the caustic property 
— we may expect the quantum abacus to operate even in 
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the semiclassical level. The important point is that, when 
we read out the outcome by the simple presence/absence 
of the particle in the qubit carrier |0) and 11} (as in the 
case of final readout of Grover algorithm 13]), the pre- 
cise profile of the qubit states does not matter so much. 
Then, this robustness of the qubit operation will further 
be improved in the semiclassical regime. 

The model considered here admits a further exten- 
sion by the addition of an extra repulsive 1/x 2 potential, 
which keeps the structure of equi-spaced energy levels in- 
tact, allowing for a quantum caustics to occur and yet 
renders the wavefunctions localized away from the gate. 
This will be useful for making the system more robust, 
especially in the instantaneous switching process between 
an open and closed gates, because the disturbance caused 
by the switching is thought to occur in the area around 
the gate. 

Once a qubit is realized, we can construct a multiple- 
qubit system by bringing many copies of our model sys- 
tem and string them together with a certain communi- 
cation mechanism to allow the multiple-qubit operation. 
The full treatment of this is outside the scope of the 
present work. We simply mention one example of two- 
qubit Control-NOT operation realized by connecting two 
of our model systems by a trigger mechanism. Namely, 
if the particle in the first system is found in x > side, 



we set a signal to travel to the second system and let the 
gate of the second system open (U = <j\) for a time step 
t, while the gate is closed (U — —I) if no signal comes 
in. The gate of the first qubit is kept closed during this 
operation. 

In mathematical terms, the two-qubit operations com- 
prise a 17(4) group, and a more direct realization may 
be to manipulate the U(4) family of point interactions 
that arise in the form of 'quantum X -junction', a graph 
of four half lines whose ends meet at a single point [lj] . 

One obvious advantage of our implementation of qubit 
operations is its visual and intuitive nature, which brings 
a pedagogical value to our model. Another advantage 
is the robustness of the qubit due to the simplicity of 
the setup, which is matched only by the spin implemen- 
tations. In contrast, most solid-state based approaches 
utilize particle states that tend to be vulnerable to 
temperature fluctuations, whose suppression can be 
costly and potentially inhibitory in constructing systems 
with a larger number of qubits. We hope that the model 
presented here offers a novel possibility for a robust and 
scalable implementation of quantum computation. 

We thank Dr. M. Iwata and Dr. Y. Kikuchi for helpful 
discussions. 
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